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Abstract 

We present in detail the calculation of the O(a^) inclusive total cross section for the 
process pp — > tih in the Standard Model, at the Tevatron center-of-mass energy y/sn = 
2 TeV. The next-to-leading order QCD corrections significantly reduce the renormalization 
and factorization scale dependence of the Born cross section. They slightly decrease 
or increase the Born cross section depending on the values of the renormalization and 
factorization scales. 



1 Introduction 



Among the most important goals of present and future colliders is the study of the electroweak 
symmetry breaking mechanism and the origin of fermion masses. If the introduction of one 
or more Higgs fields is responsible for the breaking of the electroweak symmetry, then at least 
one Higgs boson should be relatively light, and certainly in the range of energies of present 
(Tevatron) or future (LHC) hadron colliders. The present lower bounds on the Higgs boson mass 
from direct searches at LEP2 are M h > 114.1 GeV (at 95% CL) § for the Standard Model (SM) 
Higgs boson, and M h o>91.0 GeV and M A o >91.9 GeV (at 95% CL, 0.5<tan/5<2.4 excluded) 
H for the light scalar (h°) and pseudoscalar (A ) Higgs bosons of the minimal supersymmetric 
standard model (MSSM). At the same time, global SM fits to all available electroweak precision 
data indirectly point to the existence of a light Higgs boson, M h < 212 — 236 GeV 0, while 
the MSSM requires the existence of a scalar Higgs boson lighter than about 130 GeV [|]]. 
Therefore, the possibility of a Higgs boson discovery in the mass range around 115-130 GeV 
seems increasingly likely. 

In this context, the Tevatron will play a crucial role and can potentially discover a Higgs 
boson in the mass range between the present experimental lower bound and about 180 GeV 
||. The dominant Higgs production modes at the Tevatron are gluon-gluon fusion (gg — > h) 
and the associated production with a weak boson (qq — > Wh, Zh). Because of small event 
rates and large backgrounds, the Higgs boson search in these channels is extremely difficult, 
requiring the highest possible luminosity. It is therefore important to investigate all possible 
production channels, in the effort to fully exploit the range of opportunities offered by the 
available statistics. 

Recently, attention has been drawn to the possibility of detecting a Higgs signal in associ- 
ation with a pair of top-antitop quarks at the Tevatron, i.e. in pp — > tih ||. This production 
mode can play a role over most of the Higgs mass range accessible at the Tevatron. Although 
it has a small event rate, ~ 1 — 5 fb for a SM-like Higgs boson, and even lower for an MSSM 
Higgs boson, the signature (W + W~bbbb) is quite spectacular. Furthermore, at the Tevatron, 
after fully reconstructing both top quarks, the shape of the invariant mass distribution of the 
remaining bb pair is quite different for the signal and for the background. The statistics is too 
low to allow any direct measurement of the top-quark Yukawa coupling, but recent studies [|7| 
indicate that this channel can reduce the luminosity required for the discovery of a SM-like 
Higgs boson at Run II of the Tevatron by as much as 15-20%. 

The total cross section for pp — > tih has been known at tree-level for quite some time ||. 
As for any other hadronic process, next-to-leading (NLO) QCD corrections are expected to be 
important and are crucial in order to reduce the dependence of the cross section on the renor- 
malization and factorization scales. Preliminary indications of the magnitude of the NLO QCD 
corrections can be obtained in the framework of the Effective Higgs Approximation (EHA), 
where terms of order M^/y/s and Mh/m t are systematically neglected in the computation |J. 
This approximation correctly reproduces the collinear bremsstrahlung of the Higgs boson from 
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Figure 1: Feynman diagrams contributing to the lowest order process, qq — > tth. The arrows 
indicate the momentum flow. 



the heavy top quarks. However, we expect the EHA to be more reliable at the LHC center- 
of-mass energies, for which it was originally proposed, than at the Tevatron center-of-mass 
energies. We will briefly discuss the predictions of the EHA for pp — > tih in Sec. |6|. 

We also notice that QCD corrections to the associated production of a Higgs boson with 
a pair of bb quarks, which is dominated by the gg — > bbh channel, have been computed in the 
limit of large [|10|| , by resuming the leading ln(M/j/m&) terms. However this result cannot 
be applied to the tih production of a relatively light Higgs boson at the Tevatron, both because 
the ratio M h /m t is of 0(1) and does not justify the large M h limit, and because the gg channel 
is negligible for tth production at the Tevatron. 

In this paper we present in detail the calculation of the NLO inclusive total cross section for 
pp —>■ tih, cr NLO {pp —> tih), in the Standard Model, at the Tevatron center-of-mass energy. For 
pp collisions at hadronic center-of-mass energy ^fs~^=2 TeV, more than 95% of the tree- level 
cross section comes from the sub-process qq —>■ tih. Therefore, we include only the qq — > tih 
channel when computing the tree-level total cross section, and we calculate the NLO total 
cross section by adding the complete set of virtual and real 0(a s ) corrections to qq — > tih. 
The Feynman diagrams contributing to qq — > tih at lowest order are shown in Fig. [I], while 
examples of 0(a s ) virtual and real corrections are given in Figs. The main challenge in the 
calculation of the 0(a s ) virtual corrections comes from the presence of pentagon diagrams with 
several massive external and internal particles. We have calculated the corresponding pentagon 



scalar integrals as linear combinations of scalar box integrals using the method of Ref . [ ED] . The 



real corrections are computed using the phase space slicing method, in both the double [fL2j,|13| 
and single |T^-p^] cutoff approach. This is the first application of the single cutoff phase space 
slicing approach to a cross section involving more than one massive particle in the final state. 

Numerical results for our calculation of <J NLO (pp — > tih) at the hadronic center-of-mass 
energy v /s^~ = 2 TeV have been presented in |17]]. An independent calculation of the NLO total 
cross section for pp — > tih has been performed by Beenakker et at. [I8|. The numerical results 
of both calculations have been compared and they are found to be in very good agreement. The 
0(a s ) corrections to the sub-process gg — > tth are, however, crucial for determining <J NLO (pp — > 
tih) at the LHC, since pp collisions at ■J^h — 14 TeV are dominated by the gluon-gluon initial 



state. Results for the LHC are presented elsewhere |i8|,|19 
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The outline of our paper is as follows. In Sec. ^| we summarize the general structure of the 
NLO cross section, and proceed in Sees. ^ and f| to present the details of the calculation of 
both the virtual and real parts of the NLO QCD corrections. In Sec. |5| we explicitly show the 
factorization of the initial-state singularities into the quark distribution functions, and finally 
summarize our result for the NLO inclusive total cross section for pp — > tih at the Tevatron in 
Eqs. ( jS3l) and (|85|) . Numerical results for the total cross section are presented in Sec. || Explicit 
analytic results for the scalar pentagon and the infrared- singular box integrals are presented in 
Appendices |A] and ||. Appendix contains a collection of soft phase space integrals that are 
used in the calculation of the real 0(a s ) corrections to qq — > tih with the double cutoff phase 
space slicing method. Finally, in Appendix [D] we give the explicit structure of the real gluon 
emission color ordered amplitudes that are used in the calculation of the real 0(a s ) corrections 
to qq — > tih with the single cutoff phase space slicing method. 

2 General Framework 

The inclusive total cross section for pp — > tih at 0(a z s ) can be written as: 

(?nlo(pP -> tih) = J2 J dx 1 dx 2 J r f(x 1 , n)J r J(x 2 ^)a^ LO (x 1 ,x 2 , , (1) 

ij 

where J r f ,p are the NLO parton distribution functions (PDF) for parton i in a proton/antiproton, 
defined at a generic factorization scale fMf — fM, and o % ^ LO is the 0{a z s ) parton-level total cross 
section for incoming partons i and j, composed of the two channels qq, gg —>■ iih, and renor- 
malized at an arbitrary scale fi r which we also take to be fi r = fi. Throughout this paper we 
will always assume the factorization and renormalization scales to be equal, fi r —fif — fi. The 
partonic center-of-mass energy squared, s, is given in terms of the total hadronic center-of-mass 
energy squared, s H , by s = XiX 2 s H . As explained in the introduction, we consider only the 
qq — > tih channel, summed over all light quark flavors, and neglect the gg — > tih channel, since 
the gg initial state is numerically irrelevant at the Tevatron. 

We write the NLO parton-level total cross section o"^ LO (xi, x 2 , /i) as: 

°nlo( x u x 2,v) = a 2 s (/j,)^fl 1 (x 1 ,x 2 ) + ^y^/^ (zi,z 2 ,//)j (2) 

= a^ (x 1 ,x 2 ,n) + 8a% LO (xi,x 2 ,n) , 

where is the strong coupling constant renormalized at the arbitrary scale fi r = fi, a^ Q (xi, x 2 , /i) 
is the 0(a 2 s ) Born cross section, and 5a^ LO (xi,x 2 , fi) consists of the 0(a s ) corrections to the 
Born cross section, including the effects of mass factorization (see Sec. ||). 
The Born cross section to qq — > tih is given by [^1 : 
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^mj-M 2 ^ 



2 (Am 2 - M 2 )(2m 2 - M 2 ) , 8m? 



x h + 2[ - +—- - — - - + 

Xh s z sx h 



In 



x h + P" 

Xh 



where Xh = 2Eh/y/s, Eh is the Higgs boson energy in the qq center-of-mass frame, x™™ = 
2Mh/y/s, x^ ax = 1 — Ami/ s + ^hl s ^ an d we have introduced: 

n_ \ \xi-{xr n ?}[xr x -xh} Yi 2 (4) 

P \ x vaax _ Xh + Amf/s J 1 ' 

Moreover, we have defined the Yukawa coupling of the top quark to be gt = m t /v, where 
v— (Gf\/2)~ 1 I 2 is the vacuum expectation value of the SM Higgs boson, given in terms of the 
Fermi constant Gf- 

The NLO QCD contribution, Sa^ LO (xi, x 2 , //), contains both virtual and real 0(a s ) correc- 
tions to the lowest-order cross section and can be written as the sum of two terms: 

5a^ LO (x u x 2 ,fi) = J d(PS 3 )J2\A virt (ij -> tih)\ 2 + J d(PS 4 )J2\A rea i(v ^ tih + g)\ 2 

= a^(x 1 ,x 2 ,fi) + d^ sal (x 1 ,x 2 ,fJ,) , (5) 

where \ A V i rt (ij — > tih)\ 2 and \ A rea i(ij — > tt/i + ^)| 2 are respectively the squared matrix elements 
for the O(a^) zj — > tt/i and zj — > tih+g processes, and indicates that they have been averaged 
over the initial-state degrees of freedom and summed over the final-state ones. Moreover, 
d^PSz) and d(PS^) denote the integration over the corresponding three and four-particle phase 
spaces respectively. The first term in Eq. (|5]) represents the contribution of the virtual gluon 
corrections, while the second one is due to the real gluon emission. For the qq —>■ tih sub- 
process, examples of 0(a s ) virtual and real corrections are illustrated in Figs. and their 
structure is separately explained in Sees. || and f|. 

Finally, we observe that, in order to assure the renormalization scale independence of the 
total cross section at 0{a\), f l J LO ( x i, x 2 , A 4 ) i n Eq.(U) must be of the form: 

fitLofa, x 2 , fi) = fi(x 1 ,x 2 ) + fl 3 (x 1 ,x 2 )ln f^-J , (6) 
with fi(xi,x 2 ) given by: 

fl j {xi,x 2 ) = 2 \ATTb fg{xi,x 2 ) -J2 \[ dz 1 P ik (z 1 )f^ (x 1 z 1 ,x 2 ) (7) 
I k ij p 

+ [ dz 2 P kj (z 2 )f L k (x u x 2 z 2 ) 
J 



where p=(2m t + Mh) 2 /s, P%j{z) denotes the lowest-order Altarelli-Parisi splitting function |21 



of parton % into parton j, when j carries a fraction z of the momentum of parton i, (see e.g. 
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Sec. |4.1.2|) , and bo is determined by the one-loop renormalization group evolution of the strong 
coupling constant a s : 



da s {jji) 
dln(/i 2 ) 



-b a 2 s + 0(a 3 s ) 



1 /ll 



47T 



-N 



-n 



if 



(8) 



with N = 3, the number of colors, and nif — 5, the number of light flavors. The origin of 
the terms in Eq. (^) will become manifest in Sees. ^ and |], when we describe in detail the 
calculation of both virtual and real 0(a s ) corrections to qq — > tih. 
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Figure 2: 0(a s ) virtual corrections: self-energy diagrams Si' 2 ' and S2 
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3 Virtual Corrections 

The 0(a s ) virtual corrections to the tree- level qq — > tih process consist of self-energy, vertex, 
box, and pentagon diagrams which are shown in Figs. |^-|5|. We assign incoming and outgoing 
momenta according to the following notation, 

Q(Qi) + 9(?2) -> t( Pt ) + t{p' t ) + h( Ph ) , (9) 

where the momentum flow is illustrated in Figs. [||^. If we denote by Ad { the amplitude 
associated with each virtual diagram Di, the 0{a z s ) virtual amplitude squared can then be 
written as: 

SA irt | 2 -EE(^A+^ Di ) =EE 2 ^M - (10) 

i i 

5 




Figure 4: 0(a s ) virtual corrections: box diagrams B 1 , B\ ' ' and B% 
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Figure 5: virtual corrections: pentagon diagrams Pi and P2- 

ii — / p t qk 





Figure 6: C?(a s ) real corrections: examples of initial and final state real gluon emission. 



where the index i runs over the set of all virtual diagrams, and A d LO denotes the tree-level 
amplitude for qq —>■ tih calculated in d = 4 — 2e dimensions. The lowest order amplitude 
A d LO must be computed to 0(e 2 ) in order to properly account for both the singular and finite 
contributions generated by the interference of A d Q with the single and double poles present in 
the virtual amplitudes Ad v In what follows, we denote by A LO the lowest order amplitude to 
O(e ), i.e. calculated in d=4 dimensions. Also, in the following sections, the contribution of a 
given diagram or set of diagrams to a^ rt is always to be understood as the contribution of the 
corresponding term in the sum in Eq. fllPl) . 

The calculation of the virtual diagrams has been performed using dimensional regularization, 
always in d = 4— 2e dimensions. The diagrams have been evaluated using FORM |22] and Maple, 
and all tensor integrals have been reduced to linear combinations of a fundamental set of scalar 
one-loop integrals using standard techniques |23| . The scalar integrals which give rise to either 



ultraviolet (UV) or infrared (IR) singularities have been computed analytically, while finite 
scalar integrals have been evaluated using standard packages [24]. 

Self-energy and vertex diagrams contain both IR and UV divergences. The UV divergences 
are renormalized by introducing a suitable set of counterterms. Since the cross section is a 
renormalization group invariant, we only need to renormalize the wave function of the external 
fields, the top-quark mass, and the coupling constants. We discuss the renormalization of the 



UV singularities of the virtual cross section in Sec. 3.1 



Box and pentagon diagrams are ultraviolet finite, but have infrared singularities. The IR 
poles in the 0(a s ) virtual corrections are eventually canceled by analogous singularities in 
the 0(a s ) real corrections to the tree- level cross section. We discuss the structure of the IR 



8 



singularities of the virtual cross section in Sec. ^2. The structure of the IR singularities of the 
real cross section will be the subject of Sees. |4.1| and f4~2 . 

The calculation of many of the box scalar integrals and in particular of the pentagon scalar 
integrals are extremely laborious, due to the large number of massive particles present in the 
final state and in the loop. We have evaluated the necessary pentagon scalar integrals (one for 
diagram Pi and one for diagram P2), using the method of Ref. [11|, which allows the reduction of 
a scalar five-point function to a sum of five scalar four-point functions, plus terms of 0(e) which 
can be neglected. Since this is a crucial ingredient of this calculation, we will explain in detail 
in Appendix |A] how the method of Ref. JTTJ has been applied to our case. The IR-divergent 
box scalar integrals are also collected in Appendix [B]. 



3.1 UV singularities and counterterms 

The UV singularities of the 0(a^) total cross section originate from self-energy and vertex 
virtual corrections. These singularities are renormalized by introducing counterterms for the 
wave function of the external fields {5Z^\ 8Z^ ), the top-quark mass (5m t ), and the coupling 
constants (Sg t , SZ a J. If we denote by Ajjv(Si ' ) an d A-uv{Vi ' ) the UV-divergent contri- 
bution of each self-energy (Si ) or vertex diagram (V^ ) to the virtual amplitude squared 
(see Eq. ([10|) ), we can write the UV-singular part of the total virtual amplitude squared as: 



UV I 
virt I 



ElAol 2 EW, 



+ 2 



5m 



+ K (2) ) 



V 1 Juv V 2 Juv mt 



As described earlier, we denote by |-4z, | 2 the matrix element squared of the tree-level amplitude 
for qq — > tih, computed in d — 4 dimensions. We also notice that, in writing Eq. (|TT|) , we have 
included in the top-quark self-energy the top-mass counterterm, and we have used the fact that 
the Yukawa- coupling counterterm coincides with the top-mass counterterm. 

The UV-divergent contributions due to the individual diagrams are explicitly given by: 

'N 



A 



A 



A 



uv 



A 



uv 



A 



uv 



(S^ + S{ 2) 

(u, (1) + vi 2) 



a, 

2tt 

oc, 
2n' 
a. 

2tt' 



K 
K 

M t 



-N 



2N 
2 



1 



2N 



1 

£uv 



Mr, 



1 

^uv 



3 

£uv 



(12) 



1 

~2iV 



1 

£uv 
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where Af s and At are standard normalization factors defined as: 

v.= (f^)r ( i + e) , M=(^)V(i + £ ) . (is) 

Moreover we define the needed counterterms according to the following convention. For 
the external fields, we fix the wave-function renormalization constants of the external fields 
(Z^ = 1 + SZ^ 1 i = q,t) using on-shell subtraction, i.e.: 

W~ - -(£Mt-*)£) ■ 

We notice that both (JZ^ and 5Z^\ as well as some of the vertex corrections (Vi 1 ' 2 ^ and V2 ), 
have also IR singularities. In this section we limit the discussion to the UV singularities only, 
while the IR structure of these terms will be given explicitly in Sec. \i.2\ . 

We define the subtraction condition for the top-quark mass m t in such a way that m t is the 
pole mass, in which case the top-mass counterterm is given by: 

£~(£MT-Sr)GH- 

This counterterm has to be used twice: to renormalize the top-quark mass, in diagrams S{ 1] and 
and to renormalize the top-quark Yukawa coupling. As we already noted, Auv(S^ + 5i ) 
in Eq. fli~2]) already includes the top-mass counterterm. 

Finally, for the renormalization of a s we use the MS scheme, modified to decouple the top 
quark f25fl. The first nif light flavors are subtracted using the MS scheme, while the divergences 
associated with the top-quark loop are subtracted at zero momentum: 



^ = (^)(4vr)T(l + e) 



2 11 \ 1 2/1 . / fi 

^ AT \ 1 / 1 1^. / ' 



n N \ + + i n r_ 

3 3 / e uv 3\e uv \ra£ 



(16) 



such that, in this scheme, the renormalized strong coupling constant a s evolves with n// = 5 
light flavors. 



10 



It is easy to verify that the sum of all the UV-singular contributions as given in Eq. (|TTJ) is 
finite. We also notice that the left over renormalization scale dependence, due to the mismatch 
between the renormalization scale dependence of A UV (S2) and 6(Z as ), is given by: 

and corresponds exactly to the first term of Eq. (|7|), as predicted by renormalization group 
arguments. 

3.2 IR singularities 

This section describes the structure of the IR singularities originating from the 0(a s ) virtual 
corrections. The virtual IR singularities come from the following set of diagrams: vertex 
diagrams v/ 1 ' 2 ^ and V^ 1,2 "*, box diagrams B^' 2 \ box diagrams P3 1 4 , pentagon diagrams Pi 
and P2, and from the wave function renormalization of the external fields, 8Z% and 8Z^ ■ The 
IR-singular part of the total virtual amplitude squared is then of the form: 

EKfJ 2 = El Ad 2 {a ir (v« + v^) + a ir (v?> + v p) + (s4*) IR + (sz^) m 

+ A IR (p« + Pf ) + A IR (B« + pf + P X ) + A„ (pf + Pf + P,) } , (18) 

where, as before, |*4. LO | 2 denotes the matrix element squared of the tree-level amplitude for 
qq — > tih, in d = 4 dimensions. The IR-divergent contributions of the various diagrams to the 
virtual amplitude squared are given in the following: 
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where J\f s and J\f t are given in Eq. ([13]). Moreover, we have introduced the following kinematic 
invariants: 

s= s qq =2q v q 2 , s t t = 2p t -p' t , s qt = 2q v p t , s q i=2q r p' t , s qt = 2q 2 -p t , s q i = 2q T p' t , (20) 
and we have defined 



A ft - 



4m| 



1 + 



(21) 



In 



1 - & 



Substituting the explicit expression for the IR-divergent contributions given in Eq. flT9| ) into 
Eq. (H) yields: 



ZR |2 
virt I 



27 



yvirt 
T T U mH 



IR 



-IR 



(22) 



where 



- N 



1 

N 



+ 



iV 
1 

TV 



^ + ln|^]+ln|^ 



qt 



2 

-In 



5 

+ - - 



Stt 



mi 



2 (2m* + su)(3u tt 



A« - 2 In 



SqtSqt 
^qiSqt, 



(23) 



while £^y t is a finite term that derives from having factored out a common factor J\f t , and is 
given by: 



SjlR 
virt 



N 

N 



- In I -4 



1 



-In 2 — ^ 



ml 



(24) 



In Sees. |4.1.1| and [4.2.3| we will show how the IR singularities of the real cross section exactly 
cancel the IR poles of the virtual cross section (see Eqs. (|33l)- (l3~4]) and Eqs. (fT7|)- (l78|) ), as 
predicted by the Bloch-Nordsieck pH] and Kinoshita-Lee-Nauenberg [B^,E8[ theorems. 



4 Real Corrections 

The 0(a s ) corrections to qq — > tih due to real gluon emission (see Fig. |) give origin to 
IR singularities which cancel exactly the analogous singularities present in the 0(a s ) virtual 



corrections (see Sec. 3.2). These singularities can be either of soft or collinear nature and can 
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be conveniently isolated by slicing the qq — > i£/i + g phase space into different regions defined 
by suitable cutoffs, a method which goes under the general name of Phase Space Slicing (PSS). 
The dependence on the arbitrary cutoff(s) introduced in the process is not physical, and, in 
fact, cancels at the level of the total real gluon emission hadronic cross section, i.e. in a rea i, the 
real part of o NLO . This constitutes an important check of the calculation. 
We have calculated the cross section for the process 



q(qi) + 9(92) -> t{p t ) + t(p' t ) + h{p h ) + g{k) 



(25) 



using two different implementations of the PSS method which we call the two-cutoff and one- 
cutoff method respectively, depending on the number of cutoffs introduced. The two-cutoff 
implementation of the PSS method has been originally developed to study QCD corrections 
to dihadron production |Tj| and has since then been applied to a variety of processes. A nice 
review has recently appeared Jl3 to which we refer for more extensive references and details. 
The one-cutoff PSS method has been developed for massless quarks in Ref. 
to the case of massive quarks in Ref. ]16| . 



and extended 



In the next two sections we explain in detail how we have applied the PSS method to our 



case, 



using the two-cutoff implementation in Sec. |4.1| and the one-cutoff implementation in 
Sec. PO. The results for a rea i obtained using PSS with one or two cutoffs agree within the 



statistical errors. In spite of the fact that both methods are realizations of the general idea of 
phase space slicing, they have very different characteristics and finding agreement between the 
two represents an important check of our calculation. 



4.1 Phase Space Slicing method with two cutoffs 

The general implementation of the PSS method using two cutoffs proceeds in two steps. First, 
by introducing an arbitrary small soft cutoff 5 S we separate the overall integration of the 
qq — > tih + g phase space into two regions, according to whether the energy of the gluon is soft, 
i.e. Eg < 6 s s/s/2, or hard, i.e. E g > S s ^/s/2. The partonic real cross section of Eq. (|5|) can 
then be written as: 

°~real = ° 'soft + ^hard 1 (26) 

where a so f t is obtained by integrating over the soft region of the gluon phase space, and contains 
all the IR soft divergences of o qq eal . To isolate the remaining collinear divergences from a hard-, 
we further split the integration over the hard gluon phase space according to whether the gluon 
is (o'hard/coii) or is n °t {phardinon-coii) emitted within an angle 9 from the initial-state massless 
quarks such that (1 — cosO) < 5 C , for an arbitrary small collinear cutoff 5 C : 

Chard = O'hard/coii + 0~hard/ non—coll • 

(27) 

The hard non-collinear part of the real cross section, d~hard/non-coiu is finite and can be computed 
numerically, using standard Monte-Carlo techniques. In the soft and collinear regions, the 
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integration over the phase space of the emitted gluon can be performed analytically, thus 
allowing us to isolate the IR collinear divergences of a^eai- More details on the calculation of 
a so ft and &hard are given in Sec. [4.1.1| and Sec. |4.1.2j , respectively. The cross sections describing 
soft, collinear and IR-finite gluon radiation depend on the two arbitrary parameters, d~ s and S c . 
However, in the real hadronic cross section a rea i, after mass factorization, the dependence on 
these arbitrary cutoffs cancels, as will be explicitly shown in Sec. [|. 



4.1.1 Soft gluon emission 



The soft region of the qq — > tth + g phase space is defined by requiring that the energy of the 
gluon satisfies: 

V * , (28) 



Eg < 5 S 



for an arbitrary small value of the soft cutoff 5 S . In the limit when the energy of the gluon 
becomes small, i.e. in the soft limit, the matrix element squared for the real gluon emission, 
Dl^rea/I 2 ; assumes a very simple form, i.e. it factorizes into the Born matrix element squared 
times an eikonal factor $ e jfc: 

J2\ A reai(qq -> tth + g)\ 2 S -^X (ATta s )Y}A to \ 2 $ 
where the eikonal factor is given by: 

Sqt , Sq t 



eik 



(29) 



N 
4 



m. 



+ 



+ 



(pt-k) 2 (p'fk) 2 (qi-k)(p t -k) (q 2 -k)(j/fk) 



(30) 



2N 

+ *(- 



m. 



+ 



Stt 



(pt-k) 2 (p'fk) 2 (qi-k)(q 2 -k) (pt-k)(p' t -k) 

Sqt Sqt Sqt Sqt 



+ 



+ 



(q 1 -k)(p t -k) (qrk)(p' t -k) (q 2 -k)(pfk) (q 2 ■ k) (p' t ■ k) 
Moreover, in the soft region the qq — > tth + g phase space also factorizes as: 



d(PS 4 )(qq^tth + g) 



soft 



d(PS 3 )(qq^tih)d(PS g ) soft 

Md-i)u rz 
d(PS 3 )(qq - ith) , n _^Zn °( 5 ^ - E 



(31) 



(2tt)(^ 1 )2E 9 



9) i 



where d(PS g ) so f t denotes the integration over the phase space of the soft gluon. The parton 
level soft cross section can then be written as: 

a soft = (4vr« s )/i 2e J rf(PS 3 )^|Ao| 2 / d(PS g ) soft $ eik . (32) 

Since the contribution of the soft gluon is now completely factorized, we can perform the 
integration over d(PS g ) so f t in Eq. (|32|) analytically, and extract the soft poles that will have 
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to cancel X™ 1 and X™^ of Eq. (p3"l). The integration over the gluon phase space in Eq. ([32]) 
can be performed using standard techniques and we refer to Refs. [[13|,^9| for more details. For 
sake of completeness, in Appendix |C| we give explicit results for the soft integrals used in our 
calculation. 

Finally, the soft gluon contribution to a^ al can be written as follows: 



where 



X s 



-ys ris 



*'-2 



\rvirt 
~ A -2 



_ x mrt _ N 



1 

N 



+ 2 In (6, 



In 



+ -ln 2 
2 



+ 21n 2 (5 s ) -21n(<5 

r2 / c 



2, ,2 s 



1 + ln 



ml 



7l~ 

y 



In 



\m 



+ In 



SqtSqt J _ 
2 



Sill 



SqtSqt 



( ■] 



Atf 1 

+ ^ + 2 



F qt + Fi 



+ 



+ 21n(5 s ) 



In 



m? 



2 

+ ln 
1 



s 

7 2 



■In' 



2 

+2 In 



mi 



-21n 2 (5 s ) -21n(5 s 
-21nf^ 

u — - \n(— 

3 \rni 



1 + 



Sit 



(2m 2 + s t t)/3 t t 



A, 



Stt 



( SqtSgt \ 
V SqtSqt J 

2s t t 



(2m 2 t + Stt)Ptt 



Au 

" Pit 

Li 2 



2& 



1 + Af 



1, • 
+ - In 

4 



(2m 2 + s t t)Ptt 
l + Pa 



9* 



F gf +F m v-Ft 



qt 



2 In (5 S 



In 



mf 



(33) 



(34) 



while Nt is defined in Eq. (|H|), and Li 2 denotes the dilogarithm as described in Ref. |3l| . f3 t i 
and Afj are defined in Eq. (pl|) , while, for any initial parton i and final parton /, the function 
Fif can be written as: 



In 2 

,1 
+2Li 2 



1 



PfCosOif 
(3 f (l - cos% 

1-/5/ 



--In 2 
2 VI 



1 + 

0/ 



(35) 



-2Li 2 



/?/(! + cos %) 
1 — (3j cos9if 
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where cos 6if is the angle between partons i and / in the center-of-mass frame of the initial 
state partons, and 

s if 



\ 



1 



(pV 2 



1-/3/ cosOif 



(36) 



All the quantities in Eq. (^) can be expressed in terms of kinematic Al invariants, once we use 
Sif = 2pi-pf together with: 



Pt 



s - sj, + m t o _ s - s th + mf 



(37) 



where Sfh = (pf + Ph) 2 - As can be easily seen from Eqs. ( |23| ) and (pi|) , the IR poles of the 
virtual corrections are exactly canceled by the corresponding singularities in the soft gluon 
contribution. The remaining IR poles in a so f t will be canceled by the PDF counterterms as 
described in detail in Sec. B. 



4.1.2 Hard gluon emission 

The hard region of the gluon phase space is defined by requiring that the energy of the emitted 
gluon is above a given threshold. As we discussed earlier this is expressed by the condition that 



E g >5 s ^, (3> 



for an arbitrary small soft cutoff 5 S , which automatically assures that Guard does not contain 
soft singularities. However, a hard gluon can still give origin to singularities when it is emitted 
at a small angle, i.e. collinear, to a massless incoming or outgoing parton. In order to isolate 
these divergences and compute them analytically, we further divide the hard region of the 
qq — > tih + g phase space into a hard/ collinear and a hard/non- collinear region, by introducing 
a second small collinear cutoff S c . The hard/non- collinear region is defined by the condition 
that both 2 k 2 k 

J^ 1 r- > 5 C and ^ 2 ■= > 5 C (39) 

are verified. The contribution from the hard/non- collinear region, o- hard / non _ co u, is finite and 
we compute it numerically by using standard Monte Carlo integration techniques. 

In the hard/ 'collinear region, one of the conditions in Eq. ([39]) is not satisfied and the hard 
gluon is emitted collinear to one of the incoming partons. In this region, the initial-state parton 
i (i — q, q) is considered to split into a hard parton il and a collinear gluon g, i — > i'g, with 
Pi* = zpi and k — (1 — z)pi. The matrix element squared for ij — > tih + g factorizes into the 
Born matrix element squared and the Altarelli-Parisi splitting function for i i'g, i.e.: 

£i-Wy - tih + g)f C °^T (4vra s ) £ £| A LO (i'j - tih)\^-^- , (40) 

i z s ig 
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with s ig = 2pi-k. In our case: 



Piif(z) 



C F ( l -^-e{l-z) 



(41) 



is the unregulated Altarelli-Parisi splitting function for q — * q + g at lowest order, including 
terms of 0(e), and Cp = (N 2 — 1)/2N. Moreover, in the collinear limit, the qq — > tih + g phase 
space also factorizes as: 



d(PS 4 )(ij -^tth + g) 



d(PS 3 )(i'j -+ tih) { ^tZ E f - 6. 
9(cos9 ig -(l-S c )) 

d=4 -2e iXl-e) (47r) 



X 



(42) 



T(l - 2e) 16vr 2 



z ds iff [(1 — z)si S ] e ^ 



^ S 2 



where the integration range for Si g in the collinear region is given in terms of the collinear cutoff, 
and we have defined s' = 2pi> ■ pj. The integral over the collinear gluon degrees of freedom can 
then be performed separately, and this allows us to explicitly extract the collinear singularities 
of a hard . a har d/coii turns out to be of the form 



@ hard I coll 



27rr(l -2e] 



. m t 



(43) 



"l-5„ 



dz 



1 



z) 2 s' 



2z mf 



Pii'(z) o LO (i'j tth) + (i j) 



The upper limit on the z integration ensures the exclusion of the soft gluon region. As usual, 
these initial-state collinear divergences are absorbed into the parton distribution functions as 
will be described in detail in the Sec. [|]. 



4.2 Phase Space Slicing method with one cutoff 

An alternative way of isolating both soft and collinear singularities is to divide the phase space 
of the final state partons into two regions according to whether all partons can be resolved (the 
hard region) or not (the infrared, or ir, region). In the case of qq —>■ tih + g, the hard and ir 
regions are defined by whether the gluon is resolved or not. The emitted gluon is not resolved, 
and therefore considered zr, when 

s ig = 2pi ■ k < s min , with i = q,q,t,i , (44) 

for an arbitrary small cutoff s m i n . Similarly to Eq. (p6|), the partonic real cross section can be 
written as the sum of two terms: 

Klal = °ir + Vhard , (45) 
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where <7j r includes both soft and collinear singularities, while a hard is finite. Following the 
general idea of PSS, we calculate a ir analytically, while we evaluate &hard numerically, using 
standard Monte Carlo integration techniques. Both a ir and &hard depend on the cutoff s min , 
but the hadronic real cross section, a rea i, is cutoff independent, after mass factorization, as will 
be shown in Sec. |5|. 

In order to calculate we apply the formalism developed in Refs. [14-jH] as follows. 



We consider the crossed process h — > qqti + g which is obtained from qq — > tih + g by 
crossing all the initial state colored partons to the final state, while crossing the Higgs 
boson to the initial state. For a systematic extraction of the IR singularities within the 
one-cutoff method, we organize the amplitude for h — > qqti + g, A Q^* 9 , in terms of 
colored ordered amplitudes P2| . Using the color decomposition: 



rpa rpa 
1 C1C2 C3C4 



5'* 



t^c 3 e 2 jy^ciC2^c 3 C4^ , (46) 

we write A h ^ qqtt9 as the sum of four color ordered amplitudes A±, . . . , A4 as follows: 



A h^tt 9 = ig s§fqfi§ (8 ctCq T c a qCi A 1 ip t ,p' t ,q 1 ,q 2 ,k) +T^J CqCf A 2 (p t , P ' t , qi ,q 2 ,k) 



1 1 

^ctcjT^Asipup't, q u Q2, k) - —T? tC J CqC9 A 4 (pup't, Qu Q2, k) 



(47) 



where g s = ^4ira s , while (f q , fq, f t , ft) and (c q , c q , c t , Cf) denote the flavor and color indices 
of the various outgoing quarks. The amplitudes Ai(p t ,p' t , qi, q 2 , k) (for i = 1,2,3,4) 
correspond to the four possible independent color structures that arise in the h — > tiqq + g 
process, and each Ai contains terms describing the emission of the gluon from a different 
pair of external quarks. We give the explicit expressions for the Ai amplitudes in Appendix 
0. Due to this decomposition, the partonic cross section for h — > qqti + g can be written 
in a very compact form: 



with 



a 



h—*qqttg 1 2 



h^qqttg 



d(ps 5 )Y,\A 



h^qqttg |2 



4 



£{l^ 



\A 2 \ 2 + 



(48) 



(49) 



1 

iV2 



-2\A 3 + A 4 \ 2 + \A 3 \ 2 + \ Aa\ 



Using the one-cutoff PSS method and the factorization properties of both the color ordered 
amplitudes Ai and the gluon phase space in the soft /collinear limit, we extract the IR 
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singularities of a h ~' qqtt9 into o h s ^l qtt9 and d'^ qqtt9 as follows: 



&h ^u a !2IX aft** = J diPS^PS^^A^ff 9 ? , (50) 
&h ^ g corner &h ^ qi! u g = J d{ ps,) d (p Sg ) col ^\A h c - qqtla \ 2 , (51) 



where we denote by d(PS g ) so f t (d(PS g ) co u) the phase space of the gluon in the soft 
(collinear) limit, while J2\-A^ qqtt9 \ 2 {J2\A^ qqtt9 \ 2 ) represents the soft (collinear) limit of 
Eq. (f49"D . The explicit calculation of <T S( ^^f is described in detail in Sections |4.2.1| and 



[4.2.2| , respectively. The factorization of soft and collinear singularities for color ordered 
amplitudes has been discussed in the literature mainly for the leading color terms (O(N)). 
For our application of the one-cutoff PSS method, we will have to extend these results to 
the sub-leading color terms {0{1/N)). 

Finally, the IR singular contribution &i r in Eq. (^) consists of two terms: 

(Jj r (Tjy, -|- O crossing • (52) 



As described in detail in Sec. |4.2.3| , & c ir is obtained by crossing q and q to the initial state 
and h to the final state in the sum of (T sc ^f qqtt9 and cr c ^* l qqtts ', while a cross i ng corrects for the 



as 



difference between the collinear gluon radiation from initial and final state partons |L5 
will be discussed in detail in Sec. [5| As explicitly shown in Sec. OO , the IR singularities 
of a qq rt of Sec. [3]2] are exactly canceled by the corresponding singularities in <5f r . On the 



other hand, a cross i ng still contains collinear divergences that will be canceled by the PDF 
counterterms when the parton cross section is convoluted with the PDFs (see Sec. |5[). 

4.2.1 Soft gluon emission 

We first consider the case of soft singularities, when, in the limit of E g ^0 (soft limit), one or 
more Si g < s m i n (i = q, q,t,t). Using the factorization properties of the color ordered amplitudes 
Ai in the soft limit, the amplitude squared for h —> qqti + g can be written as: 

EI^-^T ^ EI<S^I 2 =(^)EI^^I 2 { fM + fM- (53) 

j-p Ifu(g) + f q q(g) - 2 (fqt{g) - fgi(g) - fqt(g) + fqt(g))] 

where, for any pair of quarks (a, b), the soft functions f a b(g) are defined as: 



/oi(ff) ,i£^_^_^, (54) 

s ag$bg S ag S b 
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and, as before (see Eq. (pOD), 

both for massless and massive quarks. A h j^ qqtt is the tree level amplitude for the process 
h — > qqtt as given by Eq. (|130|) . We note that Eq. ( |53D corresponds to the factorization 
property expressed in Eq. fl29|). Since, in the soft limit, the h qqtt + g phase space also 
factorizes, in analogy to Eq. (|3TD, we can integrate out the soft gluon degrees of freedom and 
obtain the soft gluon part of the cross section for h —>■ qqtt + g as: 



a 



h—*qqttg 
soft 



1 



d(PS/±)^2 l \A h L ^ qqtt | 2 ^ S qt + Sqt — [S U + S q q — 2 (S qt ~ S qt ~ Sq t + Sg t )} 



where, for any pair of quarks (a,b), the integrated soft functions S a b are defined as: 

9 2 S N 



S, 



ab 



d(PS g ) sof t(a,b,g)f ab (g) 



(55) 



(56) 



In the one-cutoff PSS method, the explicit form of the soft gluon phase space integral is given 
by (11: 



d(PS g ) soft (a,b,g) 



(Any aM) 



2„2 



,2 „2 



16tt 2 r(l-e 

^(Stnin S a g)6(s m i n $bg) 



SagSbgSab m b S ag m a S bg 



ds a gdsbg x 



(57) 



where 



A = < 6 - 4mX , (58) 

and the integration boundaries for s ag and Sb g vary accordingly to whether a and b are massive 
or massless quarks (see Ref. [|l6j for more details). 

The explicit form of the integrated soft functions S a b is obtained by carrying out the inte- 
gration in Eq. (p^). When a = q and b = q, i.e. when both quarks are massless, the integrated 
soft function S qq - is given by |L4| : 



2tt ) rn - e 



.5' , , 



(59) 



where, in our notation, s = s qg - is the parton center-of-mass energy (see Eq. (p0|)). On the other 
hand, when a = q,q and b = t,i, i.e. when one quark is massless and the other is massive, the 
corresponding integrated soft functions are of the form |L6] : 



S, 



ab 



2vr 
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(60) 



+ - In 2 



1 . 
+ -ln 
2 



Safe 
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Finally, when a — t and fe = t, i.e. when both quarks are massive, the corresponding integrated 
soft function S t t is given by JTB1 : 
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J s — h J a + J& 
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where we have defined: 
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and A# are defined in Eq. fl2~T|) while A# and r± are given by: 



4 - 4m 4 



t ' 



± 



to^ -1 



Finally, using Eqs. (]59[)-(|62[) , we can derive the complete form of a^^ 9 : 
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(63) 
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where 
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(64) 



4.2.2 Collinear gluon emission 

In the collinear limit when an external massless quark (i) and a hard gluon become collinear 
and cluster to form a new parton (i r ) (i + g — > z', with collinear kinematics: pj = zpj' and 
fc= (l — z)pii), the color ordered amplitudes factorize and the amplitude squared for h — > qqti+g 
can be written as: 



h^qqttg 1 2 



collinear 



ei^'t 



h^qqtt 1 2 
LO 



/r 
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><? _|_ /■?£/- 
^ •'9 
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(65) 



The collinear functions /* 9_n contain the collinear singularity and are proportional to the 
Altarelli-Parisi splitting function for ig — > i' (see Eq. (0)), i.e.: 



1-2 



(66) 



Using this definition, we can see that Eq. (p5|) is equivalent to Eq. (|40|) , although g and g, 
the massless quarks, are now considered as final state quarks. The reason why we use a more 
involved expression is because this allows us to match the collinear and soft regions of the gluon 
phase space in a very natural way, as will be explained in the following. In the same spirit, the 
lower index j of the collinear functions f] 9 ^ % keeps track of which color ordered amplitude a 
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given collinear pole comes from. Although seemingly useless at this stage, this will be crucial in 
deriving Eqs. fl6"7|) and fl6lf ), where the integration over the collinear region of the gluon phase 
space is performed in such a way to avoid to overlap with the soft gluqn phase space integration 
in Eqs. (^) and fl56|). Finally, we note that there is no fi^ 1 or f t 9 ~* in Eq. (|65| ) since the 
gluon emission from a massive quark does not give origin to collinear singularities. 

In the collinear limit, the h — > qqti + g phase space also factorizes, in complete analogy to 
Eq. (0), provided the obvious changes between initial and final state partons are taken into 
account. Therefore, we can integrate out analytically the collinear gluon degrees of freedom 
and obtain the collinear part of the partonic cross section for h — > qqti + g as: 
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(67) 



where, for any pair of quarks the integrated collinear functions CV,- are defined as: 



d(PS g ) co u(i,j,z)f j 9 ~* 1 (z) 
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The phase space of the collinear gluon can be written as: 
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d(PS g ) co u(i,j,z) 
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mm J igJ i 



(68) 
(69) 



and the integration boundaries on z are defined by the requirement that only one Si g verifies 
the condition s ig < s min . This is necessary in order to avoid overlapping with the region of phase 
space where the gluon is soft (see Eq. (^)), and it is easily translated into an upper bound on 
the z integration, thanks to the structure of Eqs. (|55|) and (|67|). In fact, each term in Eqs. (|55|) 
and fl67|) depends on only two invariants, s«, and Sj g , and each term in a^ l q ^ tth corresponds to 
an analogous term in o"^^** 9 (except that C t i is missing since there is no collinear emission 
from t and t). Therefore, for each CV,- we only need to require that when Si g <s m i n : 



jy 



> S r 



z < 1- 



s' ■ 



Z-2 



(70) 



The lower bound on z is not constrained and the integration starts at z\ = 0. For sake of 
simplicity, and since this does not give origin to ambiguities, in the following we will denote the 
s[j invariants in Eq. ( [70|) by s^-. Finally, when the integration over the collinear gluon degrees of 
freedom is performed, one finds that the /^^(^i, z 2 ) functions in Eq. are of the form ||14|| : 



Iig^i'(Zl, Zl) 




+ 0(e 2 



(71) 
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When % = q,q and j = t, t, i.e. when one quark is massless and the other is massive, the 
integrated collinear functions CV,- are given by: 

/ 47T/i 2 V 
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(72) 



while when both i,j = q, q, i.e. when both quarks are massless, 
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(73) 



Using these results, we can finally explicitly write the partonic cross section for collinear gluon 
radiation as follows: 
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4.2.3 IR-singular gluon emission: complete result for <7j, 



As already described in the beginning of Sec. [4.2| , the partonic cross section for the IR-singular 
real gluon radiation for the process qq — > tih using the one-cutoff PSS method is given by 



cr. + a, 



crossing 
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crossing • 



(76) 
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Note that crossing (y l so ^ qtt9 and d'^ u 9qtt9 only implies the interchange of the momenta of the 
quark and antiquark, since particle and antiparticle interchange under crossing. In the case of 
soft gluon emission this can be easily verified by comparing Eq. (|29"D with Eq. (|53"D , after flipping 
helicities and momenta of the crossed particles. For collinear gluon emission, the crossing is 

'"" J; - ' 1 - tlw ' difference between initial and final state collinear radiation. Using d^^ou 

), a\ r can be explicitly written as: 



in Eqs. (|63|) and 
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while Aft is defined in Eq. (|13|), and Alo is the tree-level amplitude for gg 
dimensions. 



tth in <i = 4 



As described in detail in Ref. 15 



1 ^crossing 



is given by 
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dzall(Xq^ q (z)+Xq^q(z)) 



(79) 



where X q -* q (z) (Xq^q(z)) is the unrenormalized crossing function of Ref. [ 15| , which accounts for 
the difference between collinear gluon radiation off an initial or a final state quark (antiquark): 
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5 Total cross section for pp —> tth and mass factorization 



As described in Sec. 



the observable total cross section at NLO is obtained by convoluting the 



parton cross section with the NLO quark distribution functions J-^{x,ii), thereby absorbing 
the remaining initial-state singularities of Sa^ LO into the quark distribution functions. This 
can be understood as follows. First the parton cross section is convoluted with the bare quark 
distribution functions ^^(x) and subsequently T^{x) is replaced by the renormalized quark 
distribution functions J-^(x, \x) defined in some subtraction scheme. Using the MS scheme, 
the scale-dependent NLO quark distribution functions are given in terms of T™(x) and the 



QCD NLO parton distribution function counterterms [0,[L5| as follows: 



two-cutoff PSS method 
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(82) 



where the 0(a s ) terms in the previous equations are calculated from the 0(a s ) corrections to 
the q — > qg splitting, in the PSS formalism, and P qq (z) is the Altarelli-Parisi splitting function 
of Eq. fljlD - Note that, again, we choose the factorization and renormalization scales to be 
equal. Therefore there is no explicit factorization scale dependence in Eqs. flSip and (^), and 
the only /i-dependence in ^^(x, /i) comes from a s (n). When using the two-cutoff method 
and convoluting the parton cross section with the renormalized quark distribution function of 
Eq. (PI), the IR singular counterterm of Eq. (|8lD exactly cancels the remaining IR poles of 



all t + a S0 ft and a h , 



ard/coll ■ 



In case of the one-cutoff PSS method, the IR singular counterterm 

Finally, the complete 0(a z s ) inclusive total 



of Eq. ( p2| exactly cancels the IR poles of a crossing . 
cross section for pp — ► tth in the MS factorization scheme can be written as follows: 



two-cutoff PSS method 
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Figure 7: Dependence of cr NLO (pp — > tt/i) on the arbitrary cutoff of the one-cutoff PSS method, 
s m in, at v /s^~ = 2 TeV, for M ft = 120 GeV, and ji = m t . The upper plot shows the cancellation 

and (Jhard- The lower plot shows, on an enlarged 



of the s min dependence between of , a, 



iri '-'crossmgi 



scale, the dependence of a NLO on s min , with the corresponding statistical errors. 
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We note that u NLO is finite, since, after mass factorization, both soft and collinear singularities 
have been canceled between al q irt + o' so ^ t and ahard/coii in the two-cutoff PSS method, and 
between al% t and <5f r in the one-cutoff PSS method. The last terms respectively describe the 
finite real gluon emission of Eq. (]2"7|) and Eq. (53). Note that the second term in Eqs. 
and (^), which is proportional to In (j^j, corresponds exactly to the second and third terms 
of Eq. (0), as predicted by renormalization group arguments. Before we discuss in detail the 
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Figure 8: Dependence of a NLO (pp — ► tth) on the soft cutoff S s of the two-cutoff PSS method, 
at y 7 !^" = 2 TeV, for Mh = 120 GeV, /x = m t , and <5 C = 10~ 4 . The upper plot shows the 
cancellation of the ^-dependence between a so f t + Uhard/coii and (Thard/non-coii- The lower plot 
shows, on an enlarged scale, the dependence of o^lo on S s with the corresponding statistical 
errors. 



numerical results for the NLO total cross section for pp — > tih we first demonstrate that a NLO 
does not depend on the arbitrary cutoffs of the PSS method, i.e. on s min when we use the 
one-cutoff method, and on the soft and hard/collinear cutoffs 5 S and S c when we use the two- 
cutoff method. We note that the cancellation of the cutoff dependence at the level of the 
total NLO cross section is a very delicate issue, since it involves both analytical and numerical 
contributions. It is crucial to study the behavior of ctnlo in a region where the cutoff (s) are 
small enough to justify the approximations used in the analytical calculation of the IR-divergent 
part of (J^lai, but not so small to give origin to numerical instabilities. 
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Fig. [7| is about the one-cutoff PSS method and shows the dependence of a NLO on s min . In 
the upper window we illustrate the cancellation of the s min dependence between of r , a crossing, 
and (Thard, while in the lower window we show, on a larger scale, the behavior of a NLO including 
the statistical errors from the Monte-Carlo integration. We note that a NLO also includes the 
Born cross section and the virtual contribution to the NLO cross section, which are both s m i n - 
independent, and are therefore not shown explicitly in the upper part of Fig. [7[ Clearly a 
plateau is reached in the region O.lGeV 2 < s min < 100 GeV 2 . 
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Figure 9: Dependence of a NLO (pp — > tth) on the collinear cutoff 5 C of the two-cutoff PSS 
method, at v /sj7= 2 TeV, for M h = 120 GeV, /x = m t , and 6 S = 5 x 10~ 4 . The upper plot 
shows the cancellation of the 8 C - dependence between a so f t + <7hard/coii and Uhard/non-coii- The 
lower plot shows, on an enlarged scale, the dependence of (Jnlo on 5 C with the corresponding 
statistical errors. 




Figs. H and |9] are about the two-cutoff PSS method. In Fig. |8|, we show the dependence of 
cr NLO on the soft cutoff, S s , for a fixed value of the hard/collinear cutoff, S c = 10~ 4 . In Fig. || 
we show the dependence of a NLO on the hard/collinear cutoff, S c , for a fixed value of the soft 
cutoff, S s = 5 x 10~ 4 . In the upper window of Fig. |^(^D we illustrate the cancellation of the 6 S (8 C ) 
dependence between a so f t + (T hard / coU and Ohard/non-coiu while in the lower window we show, on 
a larger scale, a NLO with the statistical errors from the Monte-Carlo integration. As before, 
<j nlo also includes the contribution from the Born and the virtual cross sections, which are 
both cutoff-independent and are not shown explicitly in the upper parts of Figs. |8|||. For S s in 
the range 10~ 4 — 2.5 x 10~ 3 and S c in the range 10~ 5 — 10~ 3 , a clear plateau is reached and the 
NLO total cross section is independent of the technical cutoffs of the two-cutoff PSS method. 
All the results presented in the following are obtained using the two-cutoff PSS method with 
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5 S and S c in the range 10 4 — 10 3 . We have confirmed them using the one-cutoff PSS method 
with 1 < s min < 10. 

6 Numerical results 

In the following we discuss in detail our results for the NLO inclusive total cross section for 
pp — > tth, a N Lo(PP ~^ tth), as introduced in Sect. || and explicitly given by Eqs. fl33|) and (|85|) . 
Our numerical results are found using CTEQ4M parton distribution functions |33| and the 
2-loop evolution of a s (p) for the calculation of the NLO cross section, and CTEQ4L parton 
distribution functions and the 1-loop evolution of a s (p) for the calculation of the lowest order 
cross section, unless stated otherwise. The top-quark mass is taken to be m t = 174 GeV and 
a N s LO {M z ) = 0.116. 
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Figure 10: Dependence of (J LO NLO {pp — > tth) on the renormalization/factorization scale p, at 
V^ = 2 TeV, for M h = 120 GeV. 




First of all, in Fig. 



10| we show how at NLO the dependence on the arbitrary renormal- 

We use M h 



ization/factorization scale p is significantly reduced. We use Mh = 120 GeV for illustration 
purposes. We note that only for scales p of the order of 2m t + M h or bigger is the NLO result 
greater than the lowest order result at ^J~s^ = 2 TeV. 

Fig. |11| shows both the LO and the NLO total cross section for pp — > tth as a function of Mh, 
at y/s~^ = 2 TeV, for two values of the renormalization/factorization scale, p = m t and p = 2mt. 



Over the entire range of Mh accessible at the Tevatron, the NLO corrections decrease the rate 
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Figure 11: cr NLO and a LO for pp — > tt/i as functions of M ft , at 1/5^" = 2 TeV, for fi = m t and 
/i = 2m ( . 
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Figure 12: K-factor for pp —>■ tth as a function of M^, at y/I^ = 2 TeV, for fi= m t and fi = 2m t 
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Mh (GeV) 

- LV - L lL ^VJ^V J 




am ffb) 

w LO \ LKJ J 


a T n (fb) 

w LO \ LKJ J 


w NLO \ LkJ J 


120 


m t 

m t + M h /2 

2m t 
zm t + ivih 


6.8662 ± 0.0013 
5.9085 ± 0.0011 
4.8789 ± 0.0009 

^t.ZO^to ZSZ U.UUUo 


5.2843 ± 0.0008 
4.5846 ± 0.0007 
3.8252 ± 0.0006 

O.OOUU ZSZ U.UUUO 


4.863 ± 0.029 

4.847 ± 0.024 

4.691 ± 0.020 
a ki i _)_ n ni 7 

41.011 Z1Z U.U1 ( 


150 


m t 

m t + M h /2 

2m t 
2m, + M h 


3.4040 ± 0.0006 
2.8289 ± 0.0005 
2.4007 ± 0.0004 
2.0282 ± 0.0004 


2.5811 ± 0.0005 
2.1668 ± 0.0004 
1.8553 ± 0.0004 
1.5813 ± 0.0003 


2.355 ± 0.013 
2.315 ± 0.011 
2.253 ± 0.010 
2.147 ± 0.008 


180 


m t 

m t + M h /2 

2m t 
2m t + M h 


1.7605 ± 0.0003 
1.4142 ± 0.0003 
1.2326 ± 0.0002 
1.0096 ± 0.0002 


1.3153 ± 0.0002 
1.0693 ± 0.0002 
0.9390 ± 0.0001 
0.7773 ± 0.0001 


1.160 ± 0.007 
1.158 ± 0.005 
1.132 ± 0.004 
1.069 ± 0.004 



Table 1: Values of both a LO (calculated with LO a s (/x) and LO PDFs), a LO (calculated with 
NLO a a (n) and NLO PDFs), and 

®nlo f° r different values of Mh and for different renormal- 

ization/f actor ization scales [i. 



for renormalization/factorization scales fi < 2m t + Mh. The reduction is much less dramatic at 
fi = 2m t than at fi = m t , as can be seen from both Fig. |TU| and Fig. |TT[ An illustrative sample of 



results is also given in Table 1. The error we quote on our values is the statistical error of the 
numerical integration involved in evaluating the total cross section. We estimate the remaining 
theoretical uncertainty on the NLO results to be of the order of 12%. This is mainly due to: 
the left over //-dependence (about 8%), the dependence on the PDFs (about 6%), and the error 
on m t (about 7%) which particularly plays a role in the Yukawa coupling. 

The corresponding K-factor, i.e. the ratio of the NLO cross section to the LO one, 

K = ^m, (86) 



is shown in Fig. [12]. For scales /x between /i = m t and /i = 2m t , the K-factor varies roughly 
between 7^ = 0.70 and f^ = 0.95, when Mh varies in the range between 100 and 200 GeV. For 
scales of the order of /i = 2m t + Mh the K-factor is of order one and becomes larger than one 
for higher scales. Given the strong scale dependence of the LO cross section, the K-factor also 
shows a significant /x-dependence and therefore is an equally unreliable prediction. Moreover it 
is important to remember that the K-factor depends on how the LO cross section is calculated. 
We choose to calculate the LO cross section using both LO a s (/i) and LO PDFs, denoted by 
a LO in Table 1. An equally valid approach could be to evaluate the LO cross section using 
NLO a s (/i) and NLO PDFs, denoted by a LO in Table 1, in which case the K-factor would 
just represent the impact of the 0(a s ) corrections that do not originate from the running 
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of and the PDFs. Since a LO > a LO , the K-factor obtained using a LO is smaller than 

the one obtained using a LO , and it is important to match the right K-factor to the right o LO 
or a LQ . Therefore we would like to stress once more that we only discuss the K-factor as a 
qualitative indication of the impact of 0(a s ) QCD corrections, for different processes or when 
using different approaches. The physical meaningful quantity is the NLO cross section, not the 
K-factor. 

For comparison, we have estimated the K-factor also in the EHA || , and we obtain K ~ 
0.6 — 0.7, for Higgs boson masses up to 150 GeV and renormalization/factorization scales in the 
range between fi = mt and fi = 2m t + M^. As anticipated, we do not expect the EHA to give a 
quantitatively good approximation of the full pp — > tih calculation at 0(a s ), since at ^fs^=2 
TeV and for a SM Higgs boson above the experimental bound, we cannot work in the limit 
Mh, mt / \fs <C 1 or Mh/mt <C 1. Still the EHA gives a remarkably good qualitative indication 
of the fact that the first order QCD corrections may lower the LO total cross section. 

It is interesting to compare our NLO result for pp — > tih with the NLO result for pp — > ti. 
Since the Higgs boson is colorless, one would naively expect the QCD corrections to both 
processes to be of roughly the same size. Defining the NLO cross section using the NLO 
evolution of and the NLO CTEQ4M PDFs, and the LO cross section using the LO 

evolution of and the LO CTEQ4L PDFs, the K-factor for it production at V /JJ7 = 2 TeV, 

with fi = m t and m t = 174 GeV, is: 

K{pp -> ti) |, 5 = 0.98 , 

K(pp->ti) \ tot = 1.05, (87) 
where the qq label indicates that only the qq initial state is included. The size of the QCD 



corrections to pp — ► ti is thus similar in magnitude to the result obtained in Fig. [12], taking 
into account that pp —>■ tih is completely dominated by the qq channel. Of course, we do not 
expect a better agreement, since in pp — > tih an additional heavy particle is produced, and 
new contributions to the virtual corrections arise. Moreover, taking the EHA as an indication, 
one could naively expect that the radiation of a Higgs boson introduces an additional negative 
contribution. We also observe that, if we now use as LO cross section the one obtained using 
NLO a a (fjL) and NLO CTEQ4M PDFs, the two K-factors in Eq. (JS7J) increase, according to the 
comment we made above, and become: 



K{pp^tt) |,5 = 1.18 , 

K(pp->ti) \ tot = 1.24, (88) 

in agreement with the literature |^4|Q- Moreover, since the NLO cross section for pp — > ti 
is further increased by the resummation of the leading and next-to-leading logarithms arising 

1 We have compared our results with Fig. 9 of Ref. j34|], and we see very good agreement with the LO and 
the NLO curves, using m t = 175 GeV and y/s^=1.8 TeV. 
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from the threshold region dynamics, the total K-factor for pp — > ti can be as high as 1.33 for 
fi = m t . With this respect, we also note that, contrary to pp — > ti, in the threshold region for 
pp — > tth there are large negative contributions, mainly from soft gluon radiation, which are 
largely compensated by large positive contributions from hard gluon radiation at larger y/s. In 
the threshold region the Coulomb term, coming from the exchange of virtual gluons between 
the t/i external legs, is important and contributes to decrease the NLO cross section, although 
it is moderated by the behavior of the three-body phase space. In the strict threshold limit, 
the Coulomb contribution to pp — > tih goes to zero, while for ti production it is constant and 
dominates the NLO cross section. 



7 Conclusion 

The NLO inclusive total cross section for the Standard Model process pp — > tth at y/s H = 2 TeV 
shows a significantly reduced scale dependence as compared to the Born result and leads to 
increased confidence in predictions based on these results. The NLO QCD corrections slightly 
decrease or increase the Born level cross section depending on the renormalization/factorization 
scales used. The NLO inclusive total cross section for Higgs boson masses in the range accessible 
at the Tevatron, 120 < M h < 180 GeV, is of the order of 1 - 5 fb. 

The contributions to the NLO cross section resulting from real gluon emission have been 
calculated in two variations of the phase space slicing method, involving one or two arbitrary 
numerical cutoff parameters, respectively. This is the first application of the one-cutoff phase 
space slicing approach, ("s m j n "), to a cross section involving more than one massive particle in 
the final state. The correspondence between the two phase space slicing approaches is made 
explicit. The virtual contributions to the NLO cross section require the calculation of both box 
and pentagon diagrams involving several massive particles and explicit results for the integrals 
have been presented in the appendices. These techniques can now be applied to other similar 
processes. 
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A Pentagon scalar integrals 



In this appendix we review the details of the calculation of the pentagon scalar integrals that 
appear in the calculation of diagrams Pi and P 2 illustrated in Fig. |5|. Using the momentum 
flow and the notation shown in Fig. ||, the pentagon scalar integral originating from diagram 
Pi (-EOpi) can be written as: 

Z7T1 4_ ^ / « ft 

where 



EOpl = /i ~ ' (2ny NiN 2 N 3 N 4 N 5 ' (89) 



Ni 


= ^ , 




N 2 


= (k + qi) 2 , 




N 3 


= (k + qi + q 2/ 




N$ 


= (k + qi + q 2 


~ Pt) 2 ~ ™t , 


N 5 


= (k + qi + q 2 


1 \2 2 

-Pt-Ph) -m t 



(90) 



We note that we included a factor /i 4_d in the definition of the (i-dimensional scalar integrals 
in order to have them in the most convenient form for the calculation of the virtual amplitude 
squared. The pentagon scalar integral originating from diagram P 2 , E0 p2 , can be obtained from 
Eqs.(|89p and ( j9"0|) by exchanging qi <-> q 2 . Therefore in the following we limit our discussion to 
E0 p i, the generalization to E0 p2 being straightforward. 



We calculate these integrals following the method introduced by the authors of Ref . [Tl[ . To 
make contact with their notation, we denote by ki the external momenta (such that kf = m 2 ), 
by Mi the internal masses, by Pi the sum of the first % external momenta, pf = J2)=i kj, by 
Pij the difference p^ —Pj-i — pf-i = fcf + k^ +1 + ■ ■ • + kj_ x (for i < j), and finally by s^- the 
invariant masses Sij — (fcj -|- kj^j . 



The topology of the generic pentagon scalar integral is illustrated in Fig. [13], which can be 
specified to our case by identifying: 

ki — > —qi (incoming q) 

k 2 — > —q 2 (incoming q) 

k 3 — ► p\ (outgoing t) (91) 

k 4 — > p h (outgoing h) 

h — > p t (outgoing t) . 



Using the standard Feynman parameterization technique, the pentagon integral in Eq.(|8 
can be written as: 

^ = -I^(^ 2 )T(3 + e) £ -\ v X)?+r • <92) 



35 




Figure 13: Topology of the pentagon scalar integral 



where the denominator T> p \{ai) is: 



V pl (di 



and the symmetric matrix Sij is given by: 
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2 j CI 2 



it; , 



(93) 



2 V 4 J 

For our particular process, the matrix SV, has the following explicit form: 
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(94) 



(95) 



Following Ref. ||11|| , -EO p i can then be written as the linear combination of five scalar box 
integrals D0 p 1 }: 



EO 



'pi 



1 5 

^E c ^° 



(i) 

pi i 



(96) 



where each DO p l scalar box integral can be obtained from the scalar pentagon integral -EO p i of 
Eq. (|H) in the limit where one of the Feynman parameters of the internal propagators goes 
to zero (i.e. DO p l is obtained when — » 0). The five box scalar integrals we need are given in 
Sees. |A.1[ - |A.5| . The coefficients q in Eq. (|9~6"D are given by: 



E^ 1 • 



(97) 



3=1 



Using Eq. (|95D we can easily obtain them in terms of m t , M/j, and the kinematic invariants Sij. 
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The final result for the pentagon scalar integral EO p i can be written as: 

X-2 , X-\ 



E ^ = 16^ 



e 2 e 



+ Xr 



(98) 



where M t is given in Eq. while X„ 2 , X_i and X are obtained using Eqs. fl96f)-fl9~5|), and 
the results in Sees. |A.1| - |A.5| . The expression for Xo is too lengthy to be given explicitly in this 
appendix, while X_ 2 and X_i have the following compact form: 

1/1 1 2 \ 
X_ 2 = — + , (99) 



2(7 V U)\T\ UJ 2 T 2 T\T 2 

1 1 

X_i = (-A CT + A W1 +A W2 -A T1 - A T2 ) + (A T2 -A n + A a 

OT\T 2 ' ~ (TT 2 LU 2 

1 ( A Tl — A T2 + A Wl ) , 



OT\<jJ\ 

where we have defined: 



o = (?i + Q2? = s , (100) 

ri = m\ - (q x - p t ) 2 = 2q 1 -p t = s qt , 

r 2 = m 2 - (q 2 - p' t ) 2 = 2q 2 -p' t = s q i , 

&i = (Pt+Ph) 2 ~ m 2 , 

^2 = (p't+Ph) 2 ~ m 2 , 



and 



A. = >„(^) , A„=ln(^) , A^ln^J , (101) 

K W1 = In ( ) , A W2 = In ( — I ] . 

\mi) \mf J 

We discuss in the following the box scalar integrals DO®, which are used in Eq. (|96l ) to calculate 
E0 P \. The analogous box scalar integrals for E0 p2 can be obtained from the DO® by exchanging 
Qi ^ Q2 m their analytic expression. 

A.l Box scalar integral DO pi 

DOpi is obtained from the pentagon in the limit a\ — ► and corresponds to the following 
integral: 



/i 4 " d 



d d k 

(2^X2X3X4X5 
d d k 1 



D *pl = I ,o-W a 7" ,r\r a 7" (102) 



(27r) d fc 2 (fc + g 2 ) 2 [(/c + g 2 - ^) 2 - m t ] [(& + 02 - p't - Ph) 2 - m 
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I] ' 



after the momentum shift k — > k — qi has been applied to the denominators listed in Eq. (]90D . 
The part of DO^} which contributes to the virtual amplitude squared is given by: 
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pi 



16vr 2 



UJ 2 T 2 
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+ x 



where Mt is given in Eq. ([13]), while the coefficients X_ 2 , X-i, and Xq are given by: 
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1 


X- 2 






2 


X- X 


= In 


X 


= #e 



5 



6 



i 2 P| 



■ 2 R 



+ 2 1nf^±^llnf^-U2 1nf^^]ln 



- 2 Li 2 



V r 2 / Vco> 2 
T2 — T1 — u 2 



T 2 



2 Li, 



W 2 

/ UJ 2 + Tl — T 2 



T"2 



^2 



2 Li, 



Vi(cU2 + Ti - 7g) 



(103) 



(104) 



(105) 



where 



2n 



ln(^lnf^U{-Li 2 (i-Ulng)ln 



m 2 t 



-ri - A+(r 2 - T X ) 



- Li 2 



T"2 



A + (r 2 - ri) + ti 



+ Li 2 



Tl 



A+(r 2 - ri) + n 



and 



A ± = -ll± 



\ 



4mf 



t 2 - n 



(A- 



(106) 



(107) 



(2) 

A. 2 Box scalar integral D0 p ± 

(2) 

D0 pl is obtained from the pentagon in the limit a 2 —> and corresponds to the following 
integral: 



DO 



(2) 
pi 



A-d 



d d k 



A-d 



(27r) d N 1 N 3 N 4 N 5 
d d k 



1 



(27i) d k 2 (k + q) 2 [(k + q- p' t ) 2 - mf] [(k + q — p't - p h ) 2 - m 2 } 



D0 p 2 2 is equal to D0 p i' , and they both coincide with DO^i in Appendix |B.l 
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,(2) 



'108) 



(3) 

A. 3 Box scalar integral DOpi 

(3) 

DOpi is obtained from the pentagon in the limit a 3 — > and corresponds to the following 
integral: 



/i 4 " d 



d d k 

{^) d NiN 2 NiN 5 
d d k 



D0$ = a A ~ d I (109) 

P 1 r" I (c_\d AT AT AT AT V / 



(2n) d k 2 (k + qi) 2 [(k + q x - p t ) 2 - m 2 }[(k + q 1 - p t - p h ) 2 - mf\ 
after the momentum shift k — > —A; — ft has been applied to the denominators in Eq. (|90D. We 



notice that this integral can be obtained from DoO* when q 2 — > qi and p' t — > p t . 



A. 4 Box scalar integral -DO^ 

-DOpi is obtained from the pentagon in the limit — ► and corresponds to the following 
integral: 



d d k 



(27r) d k 2 (k + q 2 ) 2 (k + q l + q 2 ) 2 [(k + qi + q 2 - p t ) 2 - m 



21 ' 



after the momentum shift k —k — q\ — q 2 has been applied to the denominators in Eq. 
This integral coincides with D0 b3 in Appendix [B.3| . 

(5) 

A. 5 Box scalar integral -DOpi 

(5) 

-DOpi is obtained from the pentagon in the limit — ► and corresponds to the following 
integral: 

r>n(5) _ 4-d / 1 mn 



(27r) <1 fc 2 (fc + ft) 2 + ft + ft) 2 p + ft + ft - pQ 2 - ™ 2 
This integral coincides with -DO^ in Appendix p.3| . 
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B Box scalar integrals 

B.l Box 1 : box scalar integral -DO51 

The scalar box integral DO^i can be written as: 

d d k 



D0 bl = u 4 ~ d I -^-Z- , (112) 



where 

Ni = k 2 



N 2 = (k + q) 2 , (113) 

N 3 = (k + q - p' t ) 2 - m 2 , 

iV 4 = (k + q- p' t - Vhf ~ m 2 . 



The analytical expression for this integral can be found in Ref . [35] . Since the integral is finite 



we have evaluated it in d — 4 dimensions using the FF package |24 ] . 



B.2 Box 2: box scalar integrals DO^ and DO^ 



The scalar box integral D0 b2 can be written as: 

d d k 



where 



1 v n. 

Jl) _ ,4-d 



D0$ = u 4 - rf / -ZZ- , (114) 

''- ' ' (2vr) d iV^A^ ' V ; 



iVi = k 2 , 

N 2 = (k + Pt ) 2 -m 2 , 

N 3 = (k + p t + Ph) 2 - m 2 , 

N 4 = (k-p' t ) 2 -m 2 t , (115) 

while D0%2 is obtained from DO^ by exchanging p t <-> p\. Therefore, all the following results 
for DOtf can be easily extended to DO^ ■ 

The part of DO^ which contributes to the virtual amplitude squared is of the form: 

,(1) » ,r fX-i 



mi 1 = -—Af t — i + x , (lie) 

16tt 2 V e / (m| - sthjsuPtt 
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where M t is given in Eq. while s t t = {pt + v'tf = (? - p h ) 2 > 0, (3 t t 

s t h = (Pt+Ph) 2 - The pole part X_i is: 



'1- 



4mf 



X_ x = In 



1 + 



while the finite part can be calculated using Ref. [|36 



(117) 



B.3 Box 3: box scalar integrals DO^ , D0 { ^ , D0$ , and D0$ 

The scalar box integral D0^ can be written as: 

d d k 1 



DO 



(i) 

63 



4-d 



(2vr) d N1N2N3N4 



lis) 



where 



iV 2 
iV 4 



fc 2 , 

(fc + gi) 2 , 
(A; + gi + g 2 ) 2 , 
(fc + qi + Q2 - Pt) 2 



m. 



(119) 



DQfo is obtained from DO^ by exchanging qi q 2 . On the other hand, DO^ arises from the 
box diagram where the Higgs boson is emitted from the antitop quark and corresponds to: 



N 2 
N 3 
N 4 



k 2 , 

(k + q 2 ) 2 , 

(k + qi + q 2 f , 

(k + q 1 + q 2 - p t ) 2 - m\ 



Therefore D0^ can be obtained from DOj^ by exchanging qi 



D0j$ is obtained from DO^V by exchanging qi <-> q 2 . We present here the case of DO 



l(3) 



(120) 

q 2 and p' t <-> p t . Finally, 



(i) 

63 



All 

other D0$ boxes, for 2 = 2,3,4 can be obtained following the simple pattern of substitutions 
explained above. 



The part of D0^ which contributes to the virtual amplitude squared is given by: 



DO 



(i) 

63 



16tt 2 



OT 2 ) 



X-1 



(121) 
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where jV t is defined in Eq. (0), the coefficients X_ 2 , and X are given by: 





3 




2 ' 




= ln U\ 




X 


2111 (m?) 





and a, r 2 , and to>i are defined in Eq. (|100[) . 



(122) 



C Phase space soft integrals 



In this appendix we collect the integrals which we have used in calculating the results in 
Eq. (^) starting from Eq. Q52]). For a more exhaustive treatment of the formalism used we 
refer to Refs. [|1^,|29|, from which the results in this appendix have been taken. 
We parameterize the soft gluon <i-momentum in the qq rest frame as: 



k = E g (l, . . . , sin 6\ sin 6 2} sin 6\ cos # 2 , cos^i) , 
such that the phase space of the soft gluon in d = A — 2e dimensions can be written as: 



d{PS, 



r(l-e) 7r e 



g ) soft 



T(l - 2e) (2tt) 3 Jo 
Then, all the integrals we need are of the form: 



d0isin 1_2£ 0i / d6 2 sm~' Ze 6 2 



2c , 



r d0! sin d ~ 3 6 1 I" d6 2 sin d ~ 4 2 
Jo Jo 



(a + b cos 



(A + 5 cos Oi + C sin X cos 2 ) 



In particular we need the following four cases. When A 2 7^ B 2 + C 2 , and 6 
(dropping terms of order O ((d — 4) 2 )): 



7T 



a(A + 5) Id- 4 



+ In 



(A + £?) 2 



In 2 



A 2 - B 2 - C 2 
A-VB 2 + C 2 \ 1 



A + 5 



-In" 



'A + VB 2 + C 2 
A-^[W~±ff 2 



+ 2Li 2 



B + ^B 2 + C r 
A-^fW+W 



-2Li 2 



■S - VB 2 + C 2 
A + 5 



(123) 



(124) 



(125) 



-a, we use 



(126) 
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while when b ^ —a we use: 
j(o,i) 



7T 



j(0,2) 



VB 2 + C 2 
-(d-4) 

2tt 



In 



A - VB 2 + C 2 



Li 2 



2V^ 2 + C 2 
A + V5 2 + W 



+ - In 2 



A - V5 2 + C 2 



A 2 -B 2 - C 2 
Finally, when A 2 = B 2 + C 2 , and b = —a, we have: 



i-i( d - 1) a ^ f A±^s!+£r 



jd.i) = 2tt— — f ^±£y /2_3 

aA d - 4 V 2A J 



(127) 



128) 



(129) 



D Color ordered amplitudes for /i — ► ggtt + g 

The tree-level amplitude for h — > q{qx)q{q2)t{pt)t{p't) is explicitly given by: 



(Ph-Pt-P't) 2 



7, ) T^vfa) 



(130) 



" w r>,-pi) 2 —? + fe-p.) 2 —? 

where is taken as incoming, while all the other momenta are outgoing. Using the color 
decomposition given in Eq. (46), we have rewritten A^ qqtt in terms of a leading color and a 
sub-leading color ordered amplitude. Both amplitudes are given by: 



v [Ph -Pt~ Pt) 



X 



(131) 



u (Pt) I 7*71 I7v> — + 7Z I u (Pt, 



V H 



{Ph - Pt) 2 ~ m \ {Ph - Pt) 2 - m 2 
1 -A° 



(pa - Pt - p't) 



l\2 tt,v 1 



where, for future purposes, we have introduced the A q g and A°.f tree-level partial amplitudes: 



A. 



0m 
99 



u(q l )Yv{q 2 ) 



(132) 



= 7 71 — — + — — 1^7 ) v(p t ) 



(Ph ~ P't) 2 - ml (p h - p t y 
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The 0(a s ) real corrections to the Born amplitude consist of the process h — > qqti + g, 
where the gluon can be emitted either from the external quark legs or from the internal gluon 
propagator. Therefore we can write A h ^ qqttg as follows: 



A h^m 9 = (ig s )5f q f-5f t f i \A P {T a T b ) CqCq T b tC _ + A% {T b T a ) CqCq T b tC _ 

+ At T b CqCq (T a T b ) c +AI T b CqCq (T b T a ) + A£ (if^T^T* 



(133) 



where e M (/c) is the polarization vector of the emitted gluon and we have defined by A% the 
part of the real amplitude corresponding to the emission of the gluon from i — q,q, t, i, g. More 
explicitly, the At? amplitudes are given by: 



AH 



AH 



AH 
^9 



a? 



2 m t 



fl(<b)y^-7^(<b)' 



9t- A 



2qi ■ k 

2q 2 -k 
1 



^v(q 2 



(Ph-Pt-Pt) 2 " 
\ 1 



A°' u 



(134) 



(Ph ~ Pt - p't) r "** 
1 



(Ph -Pt~ Pt) 



— (Vi l g pu ( y k,q 1 ,q 2 ) 



-A - 



2 tt,v 1 



9l-) A%u(Pt) 



r 



Jt+ ft + mt rf h - p" t + m t 



"7 



-k (p h - p' t ) 2 - mi 



P*h+ rf t + m u -$11+ rft+ ft + rrit 



(Ph ~ Pt) 2 - mj 1 (p h - ft - k) 2 - 



+ r 



mi 



■1* 



%Pt -k (Ph-Pt- k) 2 
2 m t 



mi 



) 



+ Jv 

+ lu 



(Ph-Pt) 2 - m 2 v ' ' 



(li + Q2) 2 
2p' t -k 



v (Pt) 



(Ph -p't-k) 2 - m 2 t (p h -p' t ) 2 - mf 
jh~ jt- I + m t -rf t -tf + m t A 
(Ph ~Pt- k) 2 - mj 2p' t ■ k J 



(?i + Q2) 



Mp't) 



where 



V 3 T(K gi, g 2 ) = (-2F - q p )g^ + (2q» + W + {~<f - k v )g 



hp\ 



(135) 



Using the color decomposition given in Eq. fl46|), we can also rewrite A h ^ qqtt9 as a linear 
combination of four color ordered amplitudes, as already given in Eq. P7|). By matching 
the color factors in Eq. ( |133| ) to the color factors in Eq. (f4"TD, we see that the color ordered 



amplitudes Ai(qi, q 2 , Pt, Pt, k) (for i = 1, . . . , 4) are given by |32 



A 1 (q 1 ,q 2 , Pt ,p' t ,k) = (A^ + A^-A^-e^k) 



(136) 



44 



A 2 (q 1 ,q 2 ,p uP / t ,k) = (A% + A% + Afj ■ e^k) 

A 3 {qt,q 2 ,Pt,Pt,k) = (A% + A% ) • e^k) , 
A4(qi,q2,Pt,p' t ,k) = (A? +A%) ■ e M (fc) . 
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